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SUMMARY 


A matrix method is presented for determining the longitudinal- 
stability coefficients and frequency response of an aircraft from arbi- 
trary maneuvers. The method is devised so that it can be applied to 
time-history measurements of combinations of such simple quantities as 
angle of attack, pitching velocity, load factor, elevator angle, and 
hinge moment to obtain the over-all coefficients. Although the method 
has been devised primarily for the evaluation of stability coefficients 
which are of primaiy interest in most aircraft loads and stahility 
studies, it can be used also, with a simple additional computation, to 
determine the frequency-response characteristics. The entire procedure 
can be applied or extended to other problems which can be expressed by 
linear differential equations. 


INTRODUCTION 


The longitudinal characteristics of an aircraft are often related 
by a second-order linear differential equation in which the aircraft is 
assumed to have freedom in pitch and in vertical motion; changes in 
forward velocity are so small that they can be neglected. In the evalua 
tion of tail loads, the coefficients of the differential equation and 
the elevator forcing function are generally assumed to be known and the 
response is to be determined. In the evaluation of gust problems the 
response and the coefficients are assumed to be known and the forcing 
function is to be determined. By analogy in stability and control work, 
it is desirable to determine the restoring-force and damping-force 
coefficients from known forcing functions and responses. In case the 
damping is small enoxigh to obtain the rate of decay (or logarithmic 
decrement) and period from the oscillation, the required damping and 
restoring coefficients are easily computed. Models employed in rocket- 
powered and drop tests can be and usually are so ballasted that such 



2 


NACA TN 2370 


well-defined oscillations are obtained; however, the longitudinal 
oscillations of piloted airplanes ordinarily are nearly critically damped 
and this analysis procedure cannot be applied. In any case, additional 
data and analysis are required to evaluate the control-effectiyeness 
coefficients, 

Appreciable work has been done recently in the field of determining 
the frequency-response characteristics of aircraft in flight and evalu- 
ating the stability coefficients from the frequency— response data, In 
general, the methods for determining these relationships have been to 
impose actually prescribed motions such as unit steps, triangular pulses, 
or sinusoidal motions to the elevator by means of special equipment and 
then to measure the responses, The theoretical methods for reducing such 
data are usually tailored to fit the prescribed elevator motion. Refers 
ences 1 and 2 present methods of treating input and output data by Fourier 
3-fte.lysis to determine the frequency response, Compared with the direct 
sine-wave input method of evaluating the frequency response, these 
methods require less special equipment and flight time at the expense of ’ 
additional computation. For the practical application of the Fourier 
transform method, certain restrictions are placed on the nature of the 
input and the resultant output motions: the motions must start from a 

trimmed steady-state condition and, at the end of the transient period, 
must approach either the original or the new steady or quasir-steady trim 
conditions. 

In view of the complications and limitations of existing methods of 
flight evaluation of stability coefficients and frequency response, 
development of a simple and less restricted flight test and associated 
analysis was considered desirable, A matrix method for evaluating the 
longitudinal^stability coefficients of an aircraft directly from the 
input and output time histories corresponding to arbitrary control 
motions has been derived in the present paper, The frequency response 
and some of the stability derivatives may be evaluated once these coef- 
ficients are knovyn. Although this method was derived to determine the 
second^order longitudinal response of an aircraft, it can be applied to 
other systems which can be approximated by second-order differential 
equations; extension of the method to higherr-order linear systems is also 
possible, 


A]^, A2 
b 

bt 


SYMBOLS 

combinations of aerodynamic parameters (see table I) 

wing span, feet 
tail span, feet 



chord, feet 

hinge-moment ■ coefficient 

rate of change of hinge-moment coefficient with 
elevator angle (30^/36) 

lift coefficient (L/qS) 

pitching-moment coefficient of airplane without 
horizontal tail (Mb/qS^) 

pitching— moment coefficient of isolated horizontal 
tail surface 

acceleration due to gravity, feet per second per second 
hinge moment 

airplane radius of gyration about pitching axis, feet 

empirical constant denoting ratio of damping moment 
of complete airplane to damping moment produced by 
tail 

lift, po\inds 

airplane mass, slugs (W/g) 

pitching moment of airplane, foot-pounds 

airplane load factor 

dynamic pressure, pounds per square foot 
wing area, square feet 
horizontal-tail area, square feet 
time, seconds 

true velocity, feet per second 
airplane weight, pounds 
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x-i;, • length from center of gravity of airplane to aero- 

dynamic center of tail (negative for conventional 
airplanes), feet 


Ki,K2,K3,K|^, 

K^,K^,K^,Kg, 

V O If 0 If O If .0 
% >^2 


dimensional constants occurring in equations (see 
table I) 


a 

“t 

y 

e 

6 


wing angle of attack, radians 
tail angle of attack, radians 
flight-path angle, radians 
angle of pitch (a + y) 
elevator deflection, radians 


downwash angle, radians 


"It 

0 

P 

T 

CO 


tail efficiency factor (Pt/^) 

phase angle between incremental load factor and 
elevator deflection, degrees 

mass density of air, slugs per cubic foot 

dumn^r variable of integration 

elevator angular velocity, radians per second 


The notations a and 9, a and 9, and so forth, denote single 
and double differentiations, with respect to time. ' 


a 

|a| 


bar over letter represents maximum value 

bars on sides- of symbol represent absolute value 


Matrix Notation: 



rectangular matrix 
square matrix 
column matrix 
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11^1 II 

IIAII 
IIAII ' 

■Subscripts: 

i 

0 

t 


integrating matrix (see table II) 
matrix defined by equation (2h) 
transpose of ||A|| 

denotes row elements in matrix 
denotes column elements in matrix , 
tail 
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LONGITUDINAL EQUATIONS OF MOTION 
Elevator Motion 


In this section the usual longitudinal equations of motion 
following an elevator motion are derived in such a manner as to obtain 
expressions between some of the simple combinations of variables which 
are measurable; in flight: namely, angle of attack and elevator angle,' 
pitching angular velocity and elevator angle, or load factor and elevator 
angle. The usual assiimptions of linearity, small angles, and no loss in 
airspeed during the maneuver, and no flexibility are implied. 

As in reference 3, the differential equations of motion of an air- 
plane due to a given elevator deflection may be written as (see fig. 1 for 
definitions) : 


dCi 

m^V - AaqS 

da 



A6 = 0 


( 1 ) 


dC 


m 


da 



+ 


dCLt 

dot 



a3 If . . 

V da V 


da, 

d5 


A6 


dC 




"Hi 


d6 


D. q — A6 - mK = 0 
t bt 


( 2 ) 
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By use of the definitions 


A9 = Ay + Aa 

0 = 7 + a 

• • • • • • 

9 = 7 + a 




(3) 


equations (1) and (2) are reducible to the following second-order 
differential equation giving the relation between angle of attack and 
elevator angle: 


a + K-j^a + K2 Aa » 


A6 + K^6 


Oi) 


where the K's are the constants for a given set of conditions and are 
defined in table I. The coefficient represents an effective 

aerodynamic-damping coefficient; K2 represents an effective 
aerodynamic restoring -force coefficient; and represent effec- 

tive elevator-control power coefficients. 

An alternate form of equation (h), expressing the relation between 
angle of pitch and elevator angle, may be obtained by inserting rela- 
tions of equation (3) into equation (h) and noting from equation (1) 
that 


7 = A]_ Aa + A2 A6 


7 = A^a + A26 

Pt Pt 

A7 = A]^ / Aa dt + A2 I A6 dt 

^0 Jo 


> 


( 3 ) 
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where A]_ snd A2 are combinations of aerodynamic parameters defined 

in table I, The equatipn obtained after these preceding substitutions 
are made is 


nt 

6 f K^e + K 2 A0 = A6 t / A6 dt 


( 6 ) 


where (see table I) 


Kg = Ko + K1A2 t K^Ai 


Kb “ A2K2 t A^Kg 


From the following definition for load-factor increment 


it follows that 


( 7 ) 


* W/S . ^2 

Aa = - r t An -r 7^ A6 


dCj 

dd 


Ai ^ 




if W/S • • Ap 

a = - -A 


dC 

da 


^ 6 


> 


( 8 ) 


^2 V 


^ temi in equation (8) vyas found to be small and is 
^flatted in the subsequent derivation. 



8 


NACA TN 2370 


Substituting the results from equation (8) into equation (U) yields 
another form expressing the relation between measured load-factor incre- 
ment and elevator angle as 


n + Kin + K2 An = Ky A6 + K06 (9) 

where (see table l) Ky and K0 are now different forms of the effec- 
tive control power coefficients. 


Hinge Moment 

The coefficients Ki to K0 occurring in equations (U), (6), and 

(9) are those associated with the measured elevator-motion case. The 
use of the relation 


dCu 3Cu 

Ch = ^ A6 + ^ Acu 

” 36 dat ^ 


( 10 ) 


gives the solution for A6 as 


dCv 


A6 = .t — ^Ch “ 

as’ 


aCh \ ac^. ty 


( 11 ) 


The increment in tail angle of attack to be substituted in equation (11) 
is given by 


Ao^ = 


sail - - ^ £ 2 ^ 

• da da 2 m -^ri^ 


- + -i- 

fit. 


( 12 ) 


so that 


A6 = 






— 


^ OCv 

Aa^l 

_ ^ £ s ^ \ 

da da 2 m ) 


> 


— 


\ ' 



> (13) 
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In order to shorten the subsequent derivation for the hinge-moment 
case, the term K|^6 in equation (i;) and its counterparts in equa- 
tions (6) and (9) are omitted. This effect is usually small? however, 
each individual case should be examined to see whether the term warrants 
dropping . 

A substitution of the value of A6 given by equation (13) into 
equations (U), (6), and (9) gives the following three differential 
equations for the same combination of variables with Cu and its 
integral replacing A6: . " 


a + K]_°a + K 2 ° Aa = 

(Hi) 

0 + + Kg® A0 = K3®Ch + dt 

(15) 

U 0 


n + K]_°n + Kg® An = K 5 ®Cj^ 

(16) 


INTEGRAL FORM OF EQUATIONS 


Although equations (U), (6), (9) and (lU), (15), (16) could be used 
to evaluate the effective K coefficients from flight measurements of 
Aa, 9, and An together with measurements of elevator angle, stick 
force, or hinge moment, it is seen that several differentiations of the 
measured data would be required. Inasmuch as a numerical differentiation 
process is inherently more inaccurate than the corresponding integration 
process, the, preceding equations are changed and rearranged so that 
either Aa, 0, or An, which are to be the measured values, appear as 
separate quantities on one side of the equation and the operations 
on these quantities appear on the other side. In integral form the 
rearranged equations are 

Pt nt AT nt p'z 

% j Aa dt + K 2 j J Aa dv dt - J J A6 dr dt - 



0 


A6 dt = -Aa 


(17) 
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p"b p"b n’T ^ 

Ki 0 + K2 / 0 dt - M / A6 dt - K5 / / A6 dT dt = -0 (I8) 

Jo Jo ^0 '^O 


pt pt p' 

/ An dt + K2 / / 

Jo '-'0 >^0 


t 

An dr dt - Ky / / A6 dx dt - 

Jo ^0 




Ko / A 6 dt = -An 

'0 


(19) 




pt pT 

Aa dx dt - Ko° ' / . / Cjj dx dt = -Aa 

Jo '-'0 


( 20 ) 


KiO 0 + K2° r 0 dt - Ko° r Ch dt - | T 

Jo Jo ^ Jo Jo 


Cj^ dx dt = -0 


( 21 ) 


I 




An dT dt - 



dT dt = -An 


( 22 ) 


In principle to solve any one of these equations for the K 
coefficients^ it is only necessary to tabulate the recorded values of 
the two basic variables (for example^ in equation (19) the values of 
An and A6) at a number of points t]_, t2, t^, and so forth along a 

given time history and perform th^e indicated integrations from t = 0 
up to the time of the recorded value t^. A number of simultaneous 
equations containing the unknown K* s result which are then solved. The 
number of equations can vary from a minimum, in which the number of 
ordinates is equal to the number of xinknown K*s, to the case where there 
are more equations than unknowns. When the number of ordinates equals 
the number of unknown K*s, the usual methods of solving simultaneous 
equations may be used to obtain the K*s; however, when there are more 
equations than unknowns, a least-squares method is required to reduce 
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the equations. Since the best average value of the K's is obtained 
when many points along the time history are used, a least-squares pro- 
cedure is generally preferable. 

Although the integration indicated in equations (17) to (22) can 
actually be performed graphically from the time histories, it is deemed 
better to express the equations in matrix form in order to enable a 
complete numerical solution to be made. 


MATRIX FORM OF EQUATIONS 


Since the derivation in matrix form for any one of equations (17) 
to (22) is the same as for any other equation, only equation (19), 
involving measured load factor and elevator angles, is used. In matrix 
form the system of simiiltaneous equations obtained from reading the 
time history of the load factor n against elevator angle 6 in an 
arbitrary pull-up may be written 




12 


NACA TN 2370 


In shorter form this expression may be rewritten as 



where the matrix ||a] 1 is in general a rectangular matrix; that is, for 
every time tj_, one equation or one row of the matrix ||a|| is obtained. 
The individual elements of matrix 1 |a 1| are evaluated from the known 

values of incremental load factor and incremental elevator angle. As 
mentioned previously, the integration may be performed graphically but 
in the present case, use is made of the integrating matrices derived in 
reference i;. Thus, any element in the rectangular matrix (equation (23)) 

pt pt pX 

such as / An dt or / / An dx dt may be expressed in matrix 

^0 ^0 '^0 
form as follows: 




The integrating matrix ||Ci|| as derived in reference k is given 
in table II, with a time interval At = 0.1 second. 

After the elements of the matrix A. (equations (23) and (2U)) have 
been determined either by applying the integrating matrix or by graphical 
integration, the method of least squares is applied to the solution of 
the system of simultaneous equations. In matrix notation the least- 
squares solution involves multiplication of matrix A by its transpose A’ 
so that equation (2li) becomes 
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[a'A] {Ki} = {-A* Ani} 


( 26 ) 


where the matrix [A'a] wo\ild be a U U matrix for equations (18) 
and ( 19 )., Equation (26) can now be arranged to be solved directly for 

the K's by multiplying by the inverse matrix [A'A] ^ so that finally 


Kl 


n\ 





(27) 


Alternately the system of simultaneous equations represented by 
equation ( 26 ) can be solved for the values of K by any of the well- 
known methods of solving sets of . simultaneous equations, that is, by 
eliminating the variables or by using Grout's method (reference 3 )- 
The derivation in matrix form of any of the other equations from (17) 
to (22) is similar to the plan given for equation (I9) and therefore is 
not given. 


FREQUENCY RESPONSE 


As first derived by Cornell Aeronautical Laboratory (reference 6), 
the frequency response was measured by actually subjecting the airplane 
to sinusoidal elevator motions of various frequencies by means of 
specially constructed apparatus. From these results the coefficients 
^1> forth, which are significant in control and loads work, 

could be determined provided the equation of motion was assumed. 

In the present instance since the coefficients Kq and K2 are 
determined directly from the equation of motion, the corresponding 
relations are given so that the frequency response, which is significant 
in the design of stable autopilot systems, can also be determined. 
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When a sinusoidal elevator motion has been assumed then equation (9) 
becomes 


K + K]^h + K2 An = Ky6 sin cot 


( 28 ) 


where An is the load-factor increment and co is the angular velocity 
of the elevator. Since equation (28) is a linear equation with constant 
coefficients, the steady- state solutions are of the form 


n = n sin(cot + 0) 
n = nco cos(cot + 0) f 

r\ 

ii = -nco^ sin(cot +0) 


(29) 


By a substitution of these relations into equation (28) the following 
equation is obtained: 


-nco^ sin(cot + 0) + K^nco cos(cot + 0) + K 2 n sin(cot + 0) = Ky6 sin cot (30) 


which may be rewritten as 

^(k2 - aPj sin(cot + 0) + K]_nco cos(cot + 0) = K-^6 sin cot 


(31) 


or 


B sin(cot + 0 + e) = Kyd sin cot 


(32) 


where 


B = K^6 = 


n i|(K 2 - a?y 


(Kico)^ 


(33) 




= 0 = tan ^ 


-K^co 
K 2 - co2 


(3u: 


and 
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From equation (33) the amplitude ratio of load factor to elevator angle 
is seen to be 



and the phase angle at various frequencies is given by equation (3U). 

In the present case the values of K]_, K 2 , and Ky would have 

been derived from the flight measurements and the values of co would 
be assigned. 

For the measured hinge-moments case the values of and K 2 ° 

would be used instead of and K 2 , and so forth. 

DETEroilNATION OF AERODYNAMIC DERIVATIVES 


The various K coefficients determined from the measured values 
may be termed effective coefficients and include, to some extent, effects 
of some nonlinearities, elasticity and effects of other variables which 
.are omitted in the usual analysis, In addition, as may be seen from 
table I, the K coefficients are combinations of various quantities 
involving known geometric qualities, the conditions of the problem as 
well as aerodynamic derivatives. The stability coefficients given in 
table I are expressed in a form suitable to loads work. In usual 
stability calculations, these coefficients are generally expressed in 
a simpler form where the number of aerodynamic variables are reduced 
and, as a result, the coefficients are more easily approximated. 

A total of 10 aerodynamic variables — — S ^ ^t 

da ’ da ’ da’ da^ ’ d6 ’ 

dCmt aC^ dC^ 

”d6~’ a^’ ~da'’ ^ ^appear in the definitions of the coeffi- 

"b 

cients of table I. Although all the aerodynamic derivatives cannot be 
determined directly from the four basic coefficients (namely, K 2 , 

K^, and K|^), engineering approximations of the more significant deriva- 
tives can be obtained if values are assigned to either some of the more 
accurately known derivatives or to those factors having least influence 
on .the problem. 
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The factors having the least influence on the problem are K, t)^, 

dCin. 

and the derivative — — which, respectively, allow for the contribution 

d6 

of wing-fuselage damping, tail efficiency, and moment due to tail camber 
to which average values can be assigned. A representative value of K 
is 1 . 25 . Representative values of range from about 1.2 to 0.8 with 

the higher limit applying to propeller-driven airplanes operating at low 
speed and full power and the lower limit applying at high speed with the 
propeller braking. An average value for Jets or at zero thrust for 

dCiat- 

propeller-driven airplanes is about 0.9. A representative value of 

can be obtained from existing wind-tunnel data or by using theoretical 
methods; -0.3 is an average value for tail surfaces. 

Since, as may be seen from table I, is directly proportional 

dC-i 

to an effective value of this derivative can be determined directly 

d6 

from the definition of K|^. 


In order to determine consistent values of the remaining significant 


aerodynamic derivatives 


dCL 
da ’ 


dC. 


m 


da’ 




dC 


h 


and 


d^. 

d6 ’ 


further 


da ’ da’ da^ ’ dot’ 

values must be assigned to several of the remaining derivatives. The 
derivatives chosen would naturally be those for which values could be 
obtained from other sources with the greatest degree of accuracy. 


EXAMPLES 


In order to illustrate the foregoing method as well as the consistency 
of results obtained with different sets of instrumentation, typical example 
are given using data obtained from three flights (referred to as flight 1, 
flight 2, and flight 3) of a high-speed medium Jet bomber. For. flight 1, 
the method of a computation is obtained in sufficient detail to enable a 
reader not too familiar with the mathematical details to reproduce similar 
results. Flight 1 is further divided into. case I where data for An and 
A6 /are used, and case II where data for 0 and A6 are used. Refer- 
ences 7 and 8 may be consulted for introductory discussions of least- 
squares and matrix methods. 

■ Figure 2 shows the measured time histories of velocity, altitude, 
incremental elevator displacement, incremental load factor, and incremental 
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pitching velocity obtained during a push-down pull-up maneuver. By means 
of the values from figure 2, increments in load factor and elevator angle 
at 0.1-second intervals have been tabulated in colTimns 2 and 3 of 
table III. The elements of the A matrix (equations (23) and (2U)) are 
given in columns U to 7 of table III. Each element in these columns has 
been determined by performing the indicated integrations on the results 
given in columns 2 and 3. In this instance the integrations have been 
performed by use of the previously mentioned integrating matrix derived 
in reference U. This method is particularly SToitable when automatic 
computing machines are available. 

The elements of matrix A (equation (23)) which are given in 
columns U to 7 of table III indicate that with the At spacing used, 
there are 23 equations involving the four unknown values of K. In 
order to obtain the least-squares solution of these equations, the trans- 
pose ||A|| of matrix ||A|| is required. The transpose matrix is obtained 
by interchanging . the rows and columns of matrix ||A||. 

The produce of the li-row, 23-column transpose matrix by the 23-row, 
U-column original matrix yields the U-row, ii-column matrix in the coef- 
ficients of Ki. The resulting four simxiltaneous equations are then 
solved by any of the well-known methods of solving sets of simultaneous 
equations . 

By performing the preceding operations, the following values of K 
were obtained from the data listed in table III: 


Kl 


K2 




K8 


3.311;221 


7.339706 


-119.553905 


5.819025 


In order to show how well these computed values of . K represent 
the original data, they have been reinserted into equation (19) along 
with . the measured values of A6 to, determine calculated values of An. 
The computed curve is given by the dashed line in figure 2 of the plot 
of An against t. 


The same process as was used for the relations of An and A6 was 
also applied to the relations of 9 and A6 shown in figure 2. The 
tabular material corresponding to table III is not included; the values 
of K obtained, however, were as follows: 


K. 


Ko 


K 


5 


K, 


3.13167 


8.1i123 


- 7.6212 


-12.1967 
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These values of K ^when reinserted into equation (18) resulted in the 
computed curve of 9 given by the dashed curve of figure 2. 

In addition to the preceding computations', several push-down pull- 
up maneuvers, made vinder similar conditions of altitude, weight, and 
center-of-gravity positions, were analyzed to obtain the variation of 
several of the computed K's with Mach number. In this analysis only, 
the measxirements of An and A6 were used. The results obtained for 
three Mach numbers are shown in figure 3- The short parts of the curves 
shown are the expected variations in the K's. Table I shows that 

should vary linearly with speed and the other values of K should vaiy 
parabolically . The curves shown are merely guides adjusted to as to 
pass through zero and the value of K at the 0.U5 Mach number point. 

■ The values of and K 2 shown in figiire 3 were also inserted 

into equations (3ii) and (35) to determine the corresponding curves of 
frequency response. The results are given in figure It. 

In addition the values of K 2 , and frequency response for 

case I have been computed by using the definitions of table I and aero- 
dynamic derivatives obtained from wind-tunnel tests. These results are 
also shown in figures 3 and it. The aerodynamic derivatives were listed 
in an unpublished report by the North American Aviation, Inc . , and were 
obtained in the Southern California Cooperative Wind Tunnel. 


DISCUSSION 


If only the frequency response is desired, it can be determined 
without recourse to the equations of motion; however, if the stability^ 
coefficients are desired, it will be necessary to use the equations of 
motion as has been done in the present paper. For either case several 
mathematical methods are available (references 1, 2, and 6) to obtain 
these required quantities and all methods, if carried far enough, should 
yield similar results. Thus the present method is basically no more 
accurate than any other method; however, it has the advantage of simple 
instrumentation and experimental procedures but may require more extensive 
computation. 

As with other methods where linearity is a basic assmption most 
consistent results are to be expected when the maneuvers are confined to 
the angle-of -attack region where linearity exists. In order for the 
outlined mathematical procedures to succeed, the maneuvers should cover 
as much of the linear range as possible in a short period of time and 
the portion of the maneuver considered should be confined to that portion 



NACA TN 2370 


19 


where the integrals are increasing. This practice insures that the 
elements of the original matrix A are all different and that the 
subsequent least-squares matrix [a'^ is not ill-behaved. A point 
worth noting in connection with the use of the equations is that zero 
time is assumed as being at the start of the maneuver when the airplane 
is in steady flight. Since the present method is not restricted by the 
final condition, it offers the possibility of performing an analysis 
on fragments of curves with the result that any variations in the con- 
stants may be determined. In such an analysis two possibilities occur; 
(1) where the fragments considered start from a fixed initial condition 
and become successively longer, and (2) where the fragments are taken 
as consecutive. In the first case, the present method may be applied 
without any modification; in the second case, the equations must be 
altered to introduce the initial conditions for each fragment. These 
possibilities have not, however, been explored. ^ 

In the derivation given herein, lag in downwash has been included 
(see equation (2)) but xmsteady lift effects have not. References 9 
and 10 show that for the present purposes the inaccuracy of omitting 
unsteady flow effects, except downwash lag, is probably no greater than 
the inaccviracies in the original assumptions or of the experimental data. 

Other terms and other combinations of measurements might have been 
included in the derivations given, for instance, the equations are 
readily adapted to measurement of tail load and either airplane load 
factor, airplane angle of attack, or pitching angular velocity. Addi- 
tional terms may have been included to account for flexibility. Also it 
is possible, as for example in the case of the hinge-moment relations, 
to include additional terms to account for elevator moment-of-inertia 
effect, rate of elevator motion, and so forth in order to make the 
methods more inclusive. The inclusion of these further terms, however, 
generally requires additional K's to be evaluated and would only be 
justified when the assumptions implied in the basic equations of motion 
can be more closely approached and when the accuracy of measurements is 
high. Although the method has been applied herein to second-order 
differential equations, it may be extended to higher-order equations, 
with the limitation that too many integrations destroy the conditioning 
of the equations used in determining the coefficients (equation (26)) 
and make the equations difficult to work with. 

The results of the sample computations in which two different sets 
of instrumentation were used indicate an average difference between the 
respective K coefficients of about 10 percent. The use of a least- 
squares method permits calculation of a probable error, which is an indi- 
cation of how well the second-order system and the coefficients (computed 
on the basis of 0.1-second time intervals) fit the data. The expression 
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used in computing the probable error is 


0.67U5 




where is the main diagonal term of [a'J^ E is the difference 

between the computed and measured value of the variable , and N is the 
number of cases considered in the least-squares procedure. This 
probable error has been calculated for case I and case II and indicates 
an error of ±0.3 in Kq and ±0,B in K 2 for the computations in which 

the accelerometer measurements were used. These values are contrasted 
with probable' errors of ±0.1 and ±0.3 for the pitching-angular- velocity- 
measurements . These probable errors are associated with "the very small 
differences between the solid-line and dashed-line curves shown in 
figure 2. Greater accuracy may be obtained by increasing instrument 
accuracy, record-reading accuracy, and correcting original data for 
instrument errors. Further acciiracy in the method may always be 
attained by using smaller time intervals. 


The results shown in figure 3 for the three flights investigated 
give some idea of the scatter to be expected between runs as well as 
the variation of the coefficients K^_ and K 2 with Mach number. As 

might be expected from the definition, K]_ is seen to va^y linearly 

with Mach number with little scatter. On the other hand, the values 
of K 2 either indicate a linear variation with Mach number or a 
scatter about the expected parabolic variation. 

The computed values of and K 2 (fig. 3) obtained from the 

wind-tunnel data are in fair agreement with the flight-test values. 

For mar iy engineering purposes this agreement may be adequate and probably 
typical of what might be expected if wind-tunnel data were used at the 
design stage. 


CONCLUDING REMARKS 


A matrix method has been presented for determining the longitudinal- 
stability coefficients and frequency response of an aircraft from an 
analysis of arbitrary maneuvers in which simple instrumentation is used. 
Errors in instrument accuracy and probable errors due to the use of a 
least-squares method are briefly disciissed. Possible improvements in 
the method are discussed but as of the present, it appears improvements 
would be justified only for those cases where the basic assumptions are 
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closely approached and where instrument accuracy is high. The method 
is equally applicable to other problems which can be expressed by 
second-order differential equations. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., December l3, 1950. 
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Figure 1.- Sign conventions employed. Positive directions shown. 
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T/me^ sec 


Figure 2.- Time histories of velocity, altitude, incremental elevator 
displacement, incremental load factor (computed and measured), and' 
incremental pitching velocity (computed and measured) for fli^t 1 
at Mach number 0.40. 
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Figure 4.- Airplane frequency response. 
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Abstract 


A method is presented for obtaining the over-all 
longitudinal-stability coefficients and frequency 
response of an aircraft from an analysis of arbitrary 
maneuvers in which simple instrumentation is used. 
Although the method presented deals entirely with the 
aircraft, it is equally applicable to other problems 
which can be expressed by second-order differential 
equations. 
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